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Abstract The core-degenerate (CD) scenario has been suggested to be a possible pro¬ 
genitor model of type la supernovae (SNe la), in which a carbon-oxygen white dwarf 
(CO WD) merges with the hot CO core of a massive asymptotic giant branch (AGB) 
star during their common-envelope phase. However, the SN la birthrates for this sce¬ 
nario are still uncertain. We conducted a detailed investigation into the CD scenario 
and then gave the birthrates for this scenario using a detailed Monte Carlo binary pop¬ 
ulation synthesis approach. We found that the delay times of SNe la from this sce¬ 
nario are ^70Myrs—1400Myrs, which means that the CD scenario contributes to young 
SN la populations. The Galactic SN la birthrates for this scenario are in the range of 
^7.4x 10 _5 yr _1 —3.7xl0 _4 yr _1 , which roughly accounts for ~2—10% of all SNe la. 

This indicates that, under the assumptions made here, the CD scenario only contributes 
a small portion of all SNe la, which is not consistent with the results of Ilkov & Soker 
(2013). 
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1 INTRODUCTION 

Type la supernovae (SNe la) play an important role in astrophysics. Due to their high luminosities and 
uniformity, SNe la are considered to be good distance indicators to determine cosmological parameters 
(e.g., Schmidt et al. 1998; Riess et al. 1998; Perlmutter et al. 1999). Studies on measuring cosmological 
distance through SNe la indicate that the expansion of the Universe is accelerating, which implies the 
existence of dark energy (see also, e.g., Riess et al. 2007; Kuznetsova et al. 2008). In addition, SNe la 
are also crucial for the study of galactic chemical evolution as they are the main contributes of iron to 
their host galaxies (e.g., Greggio & Renzini 1983; Matteucci & Greggio 1986). 

There is a theoretical consensus that SNe la are thermonuclear explosions of carbon-oxygen white 
dwarfs (CO WDs) in binary systems (Hoyle & Fowler 1960; Nomoto et al. 1997). However, the precise 
nature of SNe la remains uncertain, especially concerning their progenitor models and explosion mech¬ 
anisms (Leibundgut 2000; Hillebrandt & Niemeyer 2000; Parthasarathy et al. 2007; Podsiadlowski et al. 
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2008; Bogomazov & Tutukov2011; Wang & Han 2012; Wang et al. 2013b). The most remarkable SN la 
properties are their apparent similarity to each other. Most of the discussions about possible progenitors 
of SNe la are mainly concentrated upon the Chandrasekhar mass model. When the WD increases its 
mass to the Chandrasekhar mass limit, it may explode as an SN la. 

Several SN la progenitor scenarios have been proposed over the past few decades, e.g., the single¬ 
degenerate (SD) scenario where the companion of the CO WD is a non-degenerate star (e.g., Whelan & 
Iben 1973; Nomoto et al. 1984; Hachisu et al. 1996; Li & van den Heuvel 1997; Han & Podsiadlowski 
2004; Meng et al. 2009,2011; Chen & Li 2009; Wang et al. 2009a, 2010,2014b; Ablimit et al. 2014), the 
double-degenerate (DD) scenario where SNe la arise from the merging of two CO WDs (e.g., Tutukov 
& Yongelson 1981; Webbink 1984; Iben & Tutukov 1984), the double-detonation scenario where a sub- 
Chandrasekhar mass WD accumulates a layer of He-rich material from a He donor star (e.g., Woosley & 
Weaver 1994; Livne & Arnett 1995; Wang et al. 2013a), and the WD—WD collision scenario where two 
WDs collide and immediately ignite (e.g., Raskin et al. 2009; Rosswog et al. 2009; Kushniret al. 2013). 
Each of the above scenarios is not completely consistent with observations at present. Observational 
evidence suggests that these scenarios may coexist (see the review by Howell 2011; Wang & Han 2012; 
Maoz et al. 2014). 

Early numerical simulations have shown that the deficiency of the DD scenario is its tendency to 
result in an accretion-induced collapse (AIC) and, ultimately, the formation of a neutron star (Saio & 
Nomoto 1985; Hillebrandt & Niemeyer 2000) Jj] In order to overcome the AIC of the DD scenario, 
the core-degenerate (CD) scenario has been proposed. In this scenario, a Chandrasekhar or super- 
Chandrasekhar mass WD is formed from the merger of a CO WD with the hot core of an AGB star 
during the planetary nebula phase or shortly after the termination of the common envelope (CE) phase 
(e.g.. Sparks & Stecher 1974; Livio & Riess 2003; Kashi & Soker 2011; Soker 2013a; Ilkov & Soker 
2012, 2013; see also Tout et al. 2008; Mennekens et al. 2010). Soker et al. (2014) suggested that the 
properties of SN 201 lfe (e.g., the carbon rich composition of the fast moving ejector and a compact 
exploding object) may be explained by this scenario. Recently, Briggs et al. (2015) also did some popu¬ 
lation synthesis of WD+AGB core merger, and found that the majority of the high field magnetic WDs 
are the carbon-oxygen type and merge within a common envelope. 

It was originally assumed that circumstellar material (CSM) would not be present in the merger 
of two WDs (e.g., Maguire et al. 2013). However, recent theoretical studies, which investigated the 
interaction of ejected material from the WDs with the interstellar medium (Raskin & Kasen 2013; Shen 
et al. 2013), have suggested that the detectable CSM in some SNe la could be explained by the DD 
scenario (see also Ruiter et al. 2013). Soker et al. (2013) argued that a prompt violent merger via the CD 
scenario can explain the properties of SN PTF 1 lkx, e.g., the multiple shells of CSM and the interaction 
of ejected material from SN with the CSM which started 59 d after the explosion of the SN. Note that 
Dilday et al. (2012) suggested that SN PTF 1 lkx can be explained by the SD scenario. 

Although the CD scenario has many advantages, as mentioned above, which might explain some 
properties of SN la diversity, many of the characteristics of the CD scenario obtained from different 
methods are not consistent with each other, especially the SN la birthrate. Ilkov & Soker (2013) recently 
calculated the expected number of SNe la in the CD scenario, and their results showed that the CD 
scenario can account for the birthrates of SNe la within the uncertainties of several processes. The 
estimated SN la birthrate for this scenario is higher than that of observed, based on their simulations, 
assuming certain values for the parameters in their model. The purpose of this paper is to study SN la 
birthrates and delay times for this scenario using a detailed binary population synthesis (BPS) approach. 
In Section 2, we describe the BPS methods for the CD scenario. In Section 3, we show the simulation 
results of the CD scenario by the BPS approach. Finally, we present a discussion and conclusions in 
Section 4. 


1 Even if the AIC can be avoided, the WD remnant of merger might lose about 0.5 Mq via a super wind from the giant-like 
structure and fail to fulfil the critical mass for the SN explosion (e.g., Willson 2007; Shen et al. 2011). 
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2 METHODS 


Adopting tested assumptions regarding the CD scenario (Soker 2013a; Ilkov & Soker 2012, 2013), we 
performed a series of Monte Carlo simulations via Hurley’s rapid binary evolution code (Hurley et al. 
2000, 2002). In each simulation, we have followed the evolution of 1 x 10 7 sample binaries, some of 
which could form WD + AGB binaries. The criteria for potential SN la progenitors for the CD scenario 
are as follows. (1) The total mass of the WD remnant of the primary (Mwd) and the mass of the AGB 
core (Mcore, secondary) during the final common envelop (CE) phase should be super-Chandrasekhar, 
i.e., Mwd + M cor e ^ 1.4 Mq. (2) M cole ^ Mwd, the core of the AGB star has a greater mass than the 
WD remnant of the primary star. (3) The WD and the AGB core merge in the subsequent CE phase. 

The binary formation channel for the CD scenario in this paper is similar to that described in Ilkov 
& Soker (2013). The mass of the primordial primary star is in the range of 2.0-6.5 M 0 , and the initial 
mass ratio between the secondary and the primary (M 2 /M 1 ) is in the range of 0.76—1.0. The primordial 
orbital separation of the binary system should be wide enough for primary to evolve into an AGB star 
(wider than 2300 Rq). The primary loses a lot of material by the wind before it fills its Roche-lobe, 
and results in the stable Roche-lobe overflow (RLOF) which occurs later. After the stage of RLOF, the 
binary system becomes a CO WD + main sequence (MS) star system. At this stage, the secondary is 
still a MS star which is more massive than the primordial primary. The WD + MS system continues to 
evolve, and then the MS secondary becomes an AGB star and fills its Roche-lobe. After this stage, the 
system enters a CE phase owing to the deep convective envelope of the AGB star and the large mass 
ratio. In the subsequent stage, if the CE cannot be ejected, the WD will merge with the core of the AGB 
star during the CE phase (see also Soker 2013b). 

In this article. Hurley’s rapid binary volution code was adopted in our BPS approach. In this code, 
several processes are taken into consideration in the mass transfer process via RLOF, e.g., dynamical 
mass transfer, nuclear mass transfer and thermal mass transfer etc. In addition, wind accretion is also 
taken into consideration in this code. For details see Sections 2.1 and 2.6 in Hurley et al. (2002). 

We used the standard energy equations (e.g., Webbink 1984) to calculate the output of the CE phase. 
The CE is ejected if 


2df 


GM‘ on M acc 

2ai 


GM' don M eDV ^ 

^Rdon 


( 1 ) 


where A is a structural parameter that depends on the evolutionary stage of the donor, Mdon is the mass 
of the donor, M acc is the mass of the accretor, a is the orbital separation, M env is the mass of the 
donor’s envelope, i?don is the radius of the donor, and the indices i and f denote the initial and final 
values, respectively. The right side of the equation represents the binding energy of the CE, the left side 
shows the difference between the final and initial orbital energy, and a ce is the CE ejection efficiency. 
In principle, we expect 0 < a ce < 1, but we often find that a ce exceeds 1 for the purpose of explaining 
the observed binaries (Han et al. 1995; Webbink 2008). There are two highly uncertain parameters (i.e., 
A and a ce ), and we combine a ce and A into a single free parameter ct ce A in this paper. We change the 
value of ct C e A to examine its influence on the birthrates and delay time of SNe la, and set it to be 0.01, 
0.1 and 1.0 (e.g., Wang et al. 2009b). 

The basic initial parameters for the Monte Carlo BPS simulations are as follows: 

(1) A constant star formation rate (SFR) of 5 M Q yr _1 over the past 14 Gyrs is adopted, or alterna¬ 
tively, it is modeled as a delta function for a single instantaneous starburst (a burst producing 10 1 o Mq 
in stars is assumed). 

(2) The initial mass function (IMF) which proposed by Miller & Scalo (1979) is adopted. 

(3) A constant mass-ratio distribution is taken (e.g., Goldberg & Mazeh 1994). 

(4) A Monte Carlo method is utilized to generate the primordial binary samples. We assume that all 
stars are members of binary systems and that the distribution of separations is constant in log a for wide 
binaries, where a is separation and falls off smoothly at small separation: 


f £^sep (a/a 0 ) m , a<a 0 , 

tr se p, & 0 5 tt ^ CLl ; 
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where a sep ~ 0.07, ao = 10 Rq, a\ = 5.75 x 10 6 Rq = 0.13 pc and mm 1.2 (Han et al. 1995). 


(5) A circular orbit is assumed for all binaries. The orbits of semidetached binaries are generally cir¬ 
cularized by the tidal force on a timescale which is much smaller than the nuclear timescale. Moreover, 
a binary is expected to become circularized during the RLOF. As an alternative, we also consider a 
uniform eccentricity distribution in the range [0, 1]. 


(6) A substantially revised version that presented by Tout et al. (1997) is used to treat RLOF in 
Hurley’s rapid binary evolution code, and the stability of the mass transfer is described with the radius- 
mass exponent which was defined by Webbink (1985). 


(7) Metallicities were chosen to be Z=0.02. 


eo 

O 



Fig. 1 The evolution of SN la birthrates with time from the CD scenario for a SFR = 
5-MQyr -1 ) with different values of a ce A. The key to the line-styles representing different 
a ce A is giving in the upper left corner. 
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3 RESULTS 



Fig. 2 The evolution of SN la birthrate with time from the CD scenario for a single starburst of 
10 1o Mq. In this figure, the spin-down time is not included in the delay time. The open circles 
and filled squares are taken from Maoz et al. (2011) and Totani et al. (2008), respectively. The 
key to the line-styles representing different a ce X is giving in the upper right corner. 


We performed three sets of simulations to systematically investigate Galactic birthrate of SNe la for 
the CD scenario by changing the model parameter to examine their influences on the final results. The 
Galactic SN la birthrate under the assumptions made here from the CD scenario is in the range of 
~7.4x 10 _5 yr _1 — 3.7xl0~ 4 yr _1 (see Fig. 1), which accounts for ^2—10% of the Galactic SN la 
birthrate observed (~3—4 xl0~ 3 yr _1 ; Cappellaro & Turatto 1997; Li et al. 2011). The birthrate in 
this paper is much lower than the value obtained from observations. In Fig. 1, we can see that the SN la 
birthrate of the CD scenario increased as the a ce A is decreased. This trend can be understood as follows: 
the final orbital separation should be smaller to insert the same amount of energy to the envelope, which 
will lead to more WD+AGB core mergers, and hence more SNe la. Fig. 2 shows the evolution of SN 
la birthrates with time from the CD scenario for a single starburst with a total mass of 1O 1O M0, and 
the SNe la occur with no appreciable delay after merging where the spin-down time is assumed to be 
negligible. From this figure, we can see that SNe la from CD scenario occur between 70Myrs and 
1400 Myrs after the starburst, which means that the CD scenario might explain some young SNe la with 
short delay times. In this figure, the shortest SN la delay time would be mainly decided by the lifetime 
of a MS star with 6.5 M 0 . 
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Fig. 3 Distribution of the combined masses of the Mwd + M core in WD+AGB systems that 
can produced SNe la. Only systems with Af W n ^ M core are included in this figure. The solid, 
dashed and dot-dashed lines show the cases of a ce X = 0.01, a ce \ = 0.1 and a ce A = 1.0, 
respectively. 



Fig. 4 Similar to Fig. 3, but for the distribution of the mass ratio between the WD and the 
AGB core, Afwn/A/ r „ re . in the WD+AGB systems. 
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Fig. 5 Density distribution in the initial mass plane of A%d and M core in WD+AGB systems 
that can produced SNe la, in which we set a ce A = 0.01. 


In Fig. 3, we display the mass distribution of the total mass M to tai = ATwd + M cole of the 
WD+AGB systems that can ultimately produce SNe la. From this figure, we can see that there is a 
peak of il% d + ATcore i n the vicinity of 1.4Mq. This trend can be understood by the initial mass func¬ 
tion (IMF) of stars (e.g.. Miller & Scalo 1979). Fig. 4 displays the distribution of Mwn/M cole with 
different a ce X. From this figure, we can see that there is a peak of / M core in the vicinity of 1.0, 

and almost all the values of M^d / M cole are above 0.8. This result is in disagreement with that of Ilkov 
& Soker (2013), in which the peak of tfw d/A/ cor( , is in the vicinity of 0.8 indicating larger AGB cores. 
This might be a result of the much higher mass transfer parameter adopted by Ilkov & Soker (2013), 
and consequently arrived at a much lower mass ratio Mwd/M core (for details see Section 4). 

Fig. 5 presents the density distribution in the initial mass plane of Mwd and M cole where Mwd ^ 
M core . Note that the density distribution is predominantly concentrated in the vicinity of the diagonal 
in this figure. This result is in disagreement with that of Ilkov & Soker (2013), which have a larger 
distribution area above diagonal indicating more massive AGB cores. For the same reason, this might 
be a result of the much higher mass transfer parameter adopted by Ilkov & Soker (2013). 



Fig. 6 Initial orbital period distribution of WD+AGB systems with different values of a ce A, 
in which these WD+AGB systems can produced SNe la. 
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Fig. 7 Similar to Fig. 6, but for the initial mass distribution of the WDs in WD+AGB systems. 


According to our BPS approach, we also present some properties of WD+AGB systems that can 
produced SNe la, which would be helpful to search for potential progenitor systems of SNe la. Fig. 6 
shows the distribution of the initial orbital periods of the WD+AGB systems for producing SNe la. We 
note that there is a peak at ~ 10 3 4 d. In Fig.7, we display the distribution of the initial masses of the 
CO WDs in WD+AGB systems. From this figure, we can see that the distribution is in the range of 
0.64 — 1.26 Mq, and that most of the masses are concentrated in the interval between ~ 0.66 Mq and 
1.0 Mq. In Fig. 8, we show the distribution of the initial masses of the AGB stars in WD+AGB systems. 
From this figure, we can see that the distribution range is 1.4 — 7.2 Mq, and that most of the masses are 
concentrated in the interval between ~ 2 Mq and 4.5 Mq. These properties will be helpful to constrain 
the progenitor scenario studied in this paper. 

In Fig. 9, we present the distribution of the orbital period of surviving WD+WD binaries which 
are resulted from WD+AGB systems for different a ce X with a uniform eccentricity distribution. From 
this figure, we can see that as a ce \ increases, the orbital period of the surviving WD+WD binaries also 
increases. This trend can be understood as follows: a lower a ce A causes more WD+AGB core mergers 
and less surviving WD+WD binaries. Meanwhile, a lower a ce X leads to closer surviving WD+WD 
binaries, the reason of which is the same as that of the afore mentioned birthrate trend. Moreover, the 
ratios of the number of surviving WD+WD systems in this figure to those that merged WD+AGB 
systems in the CD scenario for a ce X = 0.01, 0.1 and 1.0 are 2.07, 3.45 and 15.54 according to our 
simulations, respectively. In reality, only a small portion of these surviving close WD+WD systems 
would merge by gravitational waves, as in the DD scenario, or by tidal forces at later evolution. If we 
take into consideration only those surviving WD+WD systems which merge within a Hubble timescale 
(the evolutionary timescale plus gravitational wave timescale is less than the Hubble timescale), the 
ratios are 1.66, 0.75 and 0.0 for ct ce A = 0.01, 0.1 and 1.0, respectively. From this figure, we also see that 
there are many surviving WD+WD systems, which is emerged from CE ejections, with period shorter 
than 0.01 day for a ce A=0.01. The gravitational wave radiation will bring these surviving WD+WD 
systems to merge soon after CE ejection, and this is the another sub-channel of CD scenario (see Meng 
& Yang 2012). For observational and theoretical papers on WD+WD binaries, see, e.g., Badenes & 
Maoz (2012) and Toonen et al. (2012). 
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Fig. 8 Similar to Fig. 6, but for the initial mass distribution of the AGB stars in WD+AGB 
systems. 



Fig. 9 The final orbital period distribution for the surviving WD+WD binaries which are re¬ 
sulted from WD+AGB systems with an uniform eccentricity distribution. The solid, dashed 
and dot-dashed lines show the cases of a ce A = 0.01, a ce A = 0.1 and a ce X = 1.0, respec¬ 
tively. The ratios of the number of surviving systems in this figure to those that merged in 
the CD scenario for a ce A = 0.01, 0.1 and 1.0 are 2.07, 3.45 and 15.54, respectively. If only 
those surviving WD+WD systems which merge within a Hubble timescale are taken into 
consideration, the ratios are 1.66, 0.75 and 0.0 for a ce A = 0.01, 0.1 and 1.0, respectively. 
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Fig. 10 The distribution of mass transfer parameter 77 which is adopted in our code, in which 
we set o ce A = 0.01. The dashed line shows the distribution of 77 for all potential WD+AGB 
systems. The solid line represents the distribution of 77 for those WD+AGB systems which 
can result in SNe la. 


4 DISCUSSION AND CONCLUSIONS 

We found that the Galactic SN la birthrate of the CD scenario is in the range of ~ 7.4 x 10 -5 yr _1 —3.7 x 
10 _ 4 yr _1 , which accounts for about ~2—10% of the observed value. Meng & Yang (2012) recently 
obtained the SN la birthrate via a sub-channel by which the CD scenario can form SNe la when the 
merging process of double WDs occurs within about 10 5 yr after the CE phase, the contribution from 
this sub-channel of CD scenario to all SNe la is only about 0.1%. Even when the birthrate of this sub¬ 
channel of the CD scenario is added, the total contribution of the CD scenario is still less than 10%. 
In other words, under our assumptions the birthrate of SNe la due to the CD scenario is only a small 
fraction of the total observed SN la birthrate compared with the results of Ilkov & Soker (2013). In 
contrast, Ilkov & Soker (2013) claimed that the CD scenario plays a leading role in SN la formation, 
and that the birthrate induced by the CD scenario can match the observed birthrate of SNe la by adopted 
favorite values. 

Obviously, our results significantly different from those of Ilkov & Soker (2013). The main differ¬ 
ence between our work and that of Ilkov & Soker (2013) is the treatment of the mass transfer between 
the binary constituents. Ilkov & Soker (2013) calculated the new mass of the MS secondary (M 2 ne w) 
with the formula M2 new = M2 + r/(M 1 — Mwd). At this step, the primary (with initial mass Mi) has 
evolved through the AGB phase and turned into a WD (with mass d) and the secondary (with mass 
M 2 new) is still a MS star ( with initial mass M 2 ), but has accreted mass from the primary star. In order 
to conduct a quantitative comparison with the results of Ilkov & Soker (2013), we take into account the 
distribution of the value of the mass transfer parameter 77 established by our calculations (see Fig. 10). 
In Fig. 10, we can see that the values of the mass transfer parameter 77 are in the range of ~0.22—0.43 
for all potential WD+AGB systems, and ~0.33 —0.43 for those WD+AGB systems which can result in 
SNe la. The maximum value of 77 is under 0.45, which is only half the value of 77 taken by Ilkon & Soker 
(2013), ;/=0.8—0.9. From this figure, we can also see that the mass transfer process play a crucial role 
in forming SNe la in the CD scenario. 

As a consequence of the high mass transfer parameter, the calculations of Ilkon & Soker (2013) 
resulted in more massive AGB stars and larger cores (this issue has been discussed above when analysing 
Figs 4-5), and consequently reached a much higher birthrate of SNe la. Besides, Ilkov & Soker (2013) 
estimated the number of SNe la from the CD scenario with a simple population synthesis method, and 
calculated the birthrate based on an ideal initial parameter space. In fact, the true initial parameter space 
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(capable of generating WD+AGB systems which can result in SNe la) may be relatively small, and the 
resultant birthrate of SNe la will not be so high. 

Observations of several over-luminous SNe la imply that they come from the WD explosion where 
the WD has a mass exceeding the standard Chandrasekhar limit (e.g., Howell et al. 2006; Astier et al. 
2006; Hicken et al. 2007; Yamanaka et al. 2010; Tanaka et al. 2010; Yuan et al. 2010; Scalzo et al. 2010). 
Tout et al. (2008) claimed that the formation of massive rotating WDs with strong magnetic fields might 
be attributed to the merger of a WD with the core of an AGB star. In Fig. 3, we can see that most of 
the masses of the WD+AGB core was larger than 1.4 Mq. Therefore, the CD scenario might be in a 
position to account for the formation of over-luminous SNe la, and the CD scenario contribution to SN 
la birthrates might mainly manifest itself in the form of over-luminous SNe la. Note that the SD model 
of SNe la may also produce over-luminous SNe la if the WDs have been prevented from exploding by 
the effect of differential rotation (e.g., Yoon & Langer 2005; Chen & Li 2009; Hachisu et al. 2012; Wang 
et al. 2014a). 

Pakmor et al. (2010) proposed that some SNe la are the result of a violent merger of two equally 
massive WDs, in which the mass ratio is between 0.8 and 1.0. The main SN la forming mechanism of 
the CD scenario may be the violent prompt merger of a WD with the core of a AGB star because the 
mass range and the mass ratio of the WD and the AGB core meet the criteria which was proposed by 
Pakmor et al. (2010). For example, from Fig. 4, we can see that almost all the values of M\vd/ M cole 
are concentrated in the interval of 0.8—1.0. 

In this paper, we obtained an upper limit on the birthrate of SNe la based on the CD scenario by 
taking into account all the potential WD+AGB systems. Under our assumptions and parameters the 
total contribution from all the potential WD+AGB systems to all SNe la is no more than 10%. We note 
that some of our assumptions are not in consensus with others (Ilkov & Soker 2013). We also obtained 
the delay time distribution of SNe la arising from the CD scenario, the slope of which follows a power 
law of t _1 (e.g., Graur et al. 2011; Maoz et al. 2011), and the CD scenario may explain some young SNe 
la with short delay times. The young SNe la may play an important role in galactic chemical evolution 
(Scannapieco & Bildsten 2005; Aubourg et al. 2008; Mannucci et al. 2008; Wang et al. 2009b), as a 
large amounts of iron would be returned to the interstellar medium much earlier than previously studies. 
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